Explicit characterization of the identity configuration 
in an Abehan Sandpile Model 
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Since the work of Creutz, identifying the group identities for the Abehan Sandpile Model (ASM) on 
a given lattice is a puzzling issue: on rectangular portions of complex quasi-self-similar structures 
arise. We study the ASM on the square lattice, in different geometries, and a variant with directed 
edges. Cylinders, through their extra symmetry, allow an easy determination of the identity, which 
is a homogeneous function. The directed variant on square geometry shows a remarkable exact 
structure, asymptotically self-similar. 



INTRODUCTION 

Bak, Tang and Wiesenfeld (BTW) introduced the 
Abelian Sandpile Model (with a different name) in [1] as 
a simple model for self-organized criticality. The model 
has been soon generalized from the square lattice to ar- 
bitrary graphs [2] . It is a lattice automata, described in 
terms of local height variables Zi G N, i.e. the number of 
"particles" at site i, so that a configuration is given by 
set C = {zi}. Particles fall into the system through a 
random memoryless process, and, as soon as some height 
Zi is higher than the local threshold value Zi, it topples to 
the neighbouring sites, i.e. the value Zi decreases by the 
(integer) amount A^^, while, for each site j ^ i, increases 
by the (integer) amount —Aij. This relaxation process 
(the avalanche) is guaranteed to stop in a finite number 
of steps ii := J2j ^ij — ^ f^'' sites i, and a cer- 
tain technical condition is satisfied (the matrix A must 
be a strictly dissipative diffusion kernel). All the heights 
remain non-negative if this is the case in the starting 
configuration, and Zi > An. For future convenience we 
also define bf :— J^j^fj' ^^'^ assume that also these 
quantities be non-negative. 

A graph, in general directed and with multiple edges, 
is thus defined through the non-diagonal part of —A, 
seen as an adjacency matrix, while the non-zero values 
bf^ are regarded as (in- or out-coming) "connections to 
the border" (and the sites i with bj' or b^ non-zero are 
said to be "on the border"). If A is symmetric, we have 
an undirected graph, and automatically bf = b^ We just 
write bi for bf and 6~ when clear that they are equal. 

This is the case in the original BTW model, which is 
defined on a portion of the lattice. We have in this 
case Zi = Aii = 4 on each site, and A^ = — 1 on all pairs 
of neighbouring sites. If we are on a rectangular domain, 
say of sides and Ly, we have bi = 1 for vertices on the 
sides of the rectangle, and hi = 2 for the corners. 

In [1] it is described how, already for this case, the 
avalanches show an unpredictable dynamics, with power- 
law size distribution, which candidated the model as a toy 
description of many interesting features, such as the 1// 



noise. 

Afterwords, Deepak Dhar and collaborators, in a series 
of relevant works [2-4] , elucidated a deep structure of the 
Markov chain related to the general ASM: the set of stable 
configurations S = {^j{0, 1, . . . , — 1}, the ones where 
no toppling can occur, is divided into a set T of transient 
and a set R of recurrent ones, the latter being visited 
an infinite number of times by the chain, and forming a 
single ergodic basin. Given the natural partial ordering, 
C =i C iS Zi < z[ for all i, then i? is in a sense "higher" 
than T, more precisely 

^ (C,C") e r X i? : C >- C"; (1) 

VC e T 3C" e i? : C . (2) 

In particular, the maximally- filled configuration Cmax = 
{zi — 1} is in i?, and higher than any other stable config- 
uration. 

The set R has an underlying abelian structure, for 
which a presentation is explicitly constructed in terms 
of the matrix A, through the (heavy) study of its Smith 
normal form. Shortcuts of the construction and more ex- 
plicit analytical results are achieved in the special case 
of a rectangular Lx x Ly portion of the square lattice, 
and still stronger results are obtained for the case of 

Lx = Ly [4]. 

Call hi the operator which adds a particle at site i to 
a configuration C , and a.^ the formal operator which ap- 
plies hi , followed by a sequence of topplings which makes 
the configuration stable. Remarkably, the final configu- 
ration aiC is independent from the sequence of topplings, 
and also, applying two operators, the two configurations 
ajGiC and aiUjC coincide, so that at a formal level 
and aj do commute. 

More precisely, if acting on C consists of the fall of a 
particle in ?, hi, and the sequence of topplings ti^, . . . , ti^ 
on sites ^l,...,^k^ the univocal definition of and the 
commutation of and aj follow from the two facts: 



Zj > Zj 

Zi > Zi and Zj > Zj 



tjQiC — QitjC 
ti t J — t J (-^ . 



(3) 
(4) 



2 



Another consequence is that, instead of doing all the top- 
plings immediately, we can postpone some of them after 
the following a's, and still get the same result. Similar 
manipulations show that the relation 



n«7 



(5) 



holds when applied to an arbitrary configuration. 

These facts lead to the definition of an abelian semi- 
group operation between two configurations, as the sum 
of the local height variables Zi, followed by a relaxation 
process [5]: 



C7ec"=(n<)c"=(n<-)^- 



(6) 



For a configuration C, we define multiplication by a pos- 
itive integer: 



kC = Ci 



'C 



(7) 



The set R of recurrent configurations is special, as each 
operator has an inverse in this set, so the operation 
above, restricted to R, is raised to a group operation. Ac- 
cording to the Fundamental Theorem of Finite Abelian 
Groups, any such group must be a "discrete thorus" 
Zdi X X • • • X Zd , for some integers di > d2 > ■ ■ ■ > dg, 
and such that da+i divides da for each a = 1, . . . , g — 1. 
The values da, called elementary divisors of A, and a 
set of generators with the proper periodicities, can be 
constructed through the normal form decomposition [4]. 
The composition of whatever C = {zi} with the set R 
acts then as a translation on this thoroidal geometry. A 
further consequence is that, for any recurrent configura- 
tion C, the inverse configuration (— C) is defined, so that 
A; C is defined for fc G Z. 

Consider the product over sites i of equations (5) 



-Ay 



(8) 



On the set R, the inverses of the formal operators a,; are 
defined, so that we can simplify common factors in (8), 
recognize the expression for 6+, and get 



n- 



(9) 



so that Yl^a^' C = C is a necessary condition for C to 
be recurrent (but it is also sufficient, as no transient con- 
figuration is found twice in the same realization of the 
Markov chain), and goes under the name of identity test. 

This condition turns into an equivalent and computa- 
tionally-cheaper procedure, called burning test, of which 
a side product, in case of positive answer, is a spanning 
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Figure 1: Left: a portion of the square lattice. Right: a 
portion of the directed square lattice we considered in this 

work, with in- and out-odgcs alternated cyclically, and white 
and gray faces with arrows oriented clockwise and counter- 
clockwise. 



arborescence rooted on the vertices of the border. So, the 
burning test provides us a bijection between the two en- 
sembles, of recurrent sandpile configurations and rooted 
arborescences with roots on the border. This is in agree- 
ment with Kirchhoff Matrix- Tree theorem, which states 
that the number of such arborescences is given by det A, 
while the number of recurrent configurations is known to 
be det A as the first step of the procedure which deter- 
mines the elementary divisors da [3]. 

If, for the graph identified by A, a planar embedding 
exists, with all sites i such that 6^ > on the most ex- 
ternal face, then the planar dual of a rooted arborescence 
coming from the burning test is a spanning tree on the 
planar-dual graph. 

This is clear for the undirected case A = A^, and needs 
no more words. If the graph is directed, the arborescence 
is directed in the natural way, while the dual tree is con- 
strained in some complicated way (some combinations of 
edge-occupations are forced to fixed values). However, 
the graphical picture simplifies if, on the planar embed- 
ding, in-coming and out-coming edges are cyclically al- 
ternating (and all plaquettes have consistent clockwise or 
counter-clockwise perimeters) , and all sites with either 6+ 
or b~ positive are on the most external face, with in- and 
out-going arrows cyclically alternating (cfr. for example 
Figure 1). We will call a directed graph of this kind a 
planar alternating directed graph. 

The connection with uniform spanning trees and the 
Kirchhoff theorem explains a posteriori the arising of self- 
organized criticality, i.e. the appearence of long-range be- 
haviour with no need of tuning any parameter. Indeed, 
uniform spanning trees on regular 2-dimensional lattices 
are a c = — 2 logarithmic Conformal Field Theory (CFT), 
and have no parameter at all, being a peculiar limit — *• 
of the Potts model in Fortuin-Kasteleyn formulation [6], 
or a limit of zero curvature in the 0SP(1|2) non-linear 
(7- model [7] . If instead one considers the larger ensemble 
of spanning forests, in a parameter t counting the compo- 
nents (or describing the curvature of the 0SP(1|2) super- 
sphere), the theory in two dimensions is scale- invariant 
for three values: at t = (the spanning trees, or the 
endpoint of the ferromagnetic critical line of Potts), at 
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the infinite-temperature point t = oo, and at some non- 
universal negative t corresponding to the endpoint of the 
anti- ferromagnetic critical line of Potts, being tc = —1/4 
on the square lattice. Through the correspondence with 
the non-linear cr-model, one can deduce at a perturbative 
level the renormalization group flow, and, in particular, 
that the system is asymptotically free for t > [7, 8], so 
that i = 0+ is an ultraviolet fixed point. 

An interesting connection between the height variables 
in the ASM and the fields of the underlying CFT is per- 
formed in [3, 9-11], and especially in [11]. 



IDENTITIES 

Given the algebraic relation (9), and the semi-group op- 
eration (6), one could define the frame configuration Idf 
as the one with Zi = bj' for all i, and realize that it acts 
as an identity on recurrent configurations, Id[ (B C = C 
if C € R. Conversely, in general it does not leave un- 
changed a transient configuration, and in particular, as, 
for any relevant extensive graph, Id{ is itself transient, 
we have that Idf © Idf, and Idf © Idf © Idf, and so on, 
are all different, up to some number of repetitions k at 
which the configuration is sufficiently filled up with par- 
ticles to be recurrent. We call Id^ this configuration, and 
fc(A) the minimum number of repetitions of Idf required 
in the ASM identified by A (we name it the filling num- 
ber of A). The configuration Idr is the identity in the 
abelian group Z^j x • • • x Z^^ described above, and to- 
gether with a set of generators e^, completely identifies 
in a constructive way the group structure of the statis- 
tical ensemble. The relevance of this configuration has 
been stressed first by M. Crcutz [5], so that we shall call 
it Creutz identity. See fig. 2 for an example. 

Unfortunately, despite many efforts, it has not been 
possible to give a closed-formula recipe for this identity 
state on given large lattices, not even in the case of a 
L X L square, and the direct numerical investigation of 
these configurations at various sizes has produced pecu- 
liar puzzling pictures [12]. 

In a large-side limit, we have the formation of curvi- 
linear triangular regions of extensive size (of order L), 
filled with regular patterns, and occasionally crossed 
by straigth "defect lines", of widths of order 1, which, 
furthermore, occasionally meet at Y-shaped "scattering 
points", satisfying peculiar conservation laws [13]. 

Similar features emerge also for the filling numbers, 
e.g. on the square lattice of size L, the index kL = 
fc(A(^)) is not badly fitted, for even L, by a parabola 
kL — + o{L'^), but showing fluctuations due to 

unpredictable numbcr-thcorctical properties of L. The 
challenging sequence of these numbers, for L up to 64, is 
given in Table I. It should be noted that, conversely, odd 
sizes 2L + 1 are related to 2L through a property proven 
in [4, sec. 7]. 




Figure 2: The top-left corners of the identities Idi^^ for L = 

50 and L — 100 in the BTW model (the other quadrants are 
related by symmetry). Heights from to 3 correspond to gray 
tones from dark to light. The smaller-size identity is partially 
reproduced at the corner of the larger one, in a fashion which 
resembles the results of Theorem 1. 



L 


2 4 6 8 10 12 14 16 18 20 


22 


> 


kL 


1 4 7 13 19 27 35 46 58 71 


87 


► 




24 26 28 30 32 34 36 38 40 42 


44 




► 


103 119 138 156 180 198 226 248 276 305 


334 


► 


> 


46 48 50 52 54 56 58 60 62 64 




► 


367 397 430 464 499 538 572 615 653 699 





Table 1: Values of fei for the BTK Abehan Sandpile on square 
geometries of even size, for L = 2, . . . , 64. 



The determination of Idy for a given graph is a proce- 
dure polynomial in the size of the graph. E.g. one could 
prove that fc(A) is sub-exponential, and that relaxing 
C © C" for C and C both stable is polynomial, then 
one can produce the powers 2* Idf recursively in s, up to 
get twice the same configuration. Better algorithms exist 
however, see for example [5]. 

Still, one would like to have a better understanding of 
these identity configurations. This motivates us to the 
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study of the ASM on different regular graphs, such that 
they resemble as much as possible the original model, 
but simplify the problem in some regards so that the 
family (over L) of resulting identities can be understood 
theoretically. 

In this research, we keep in mind a few principles: 

(1) wc wish to simplify the problem in some structural 
way; (2) we want to preserve the property of the L x L 
lattice, discussed in [4, sec. 4], of reconstruction of the 
elementary divisors of A from a suitable matrix n.yy/, 
outcome of a machinery for producing an "economic" 
presentation of the group, in terms of 0{L) generators 
fli of the whole set of L^; (3) possibly, wc want to pre- 
serve planarity, and the interpretation of dual spanning 
subgraphs as spanning trees on the dual lattice, i.e. we 
want to use either planar undirected graphs, or planar 
alternating directed graphs (according to our definition 
above). 

A PROLOG: CYLINDRIC GEOMETRIES 

A first possible simplification comes from working on a 
cylindrical geometry. We call respectively periodic, open 
and closed the three natural conditions on the boundaries 
of a X Ly rectangle. For example, in the BTW Model, 
for site (i, 1), in the three cases of periodic, open and 
closed boundary conditions on the bottom horizontal side 
we would have the toppling rules 



periodic: 


ZiS 


Zi 


.1-4 


Zi±lA 


Zi±l,l + 1 




Zi 


,2 


Zi 


,2 + 1 


Zi,Ly 


Zt,Ly + 1; 


open: 


Zi 


,1 


— > Zi 


,1-4 


Zi±l,l 


Zi±\,\ + 1 




Zi 


,2 


— > Zi 


,2 + 1 






closed: 


Zi 


,1 


— > Zi 


,1-3 


Zi±l,l 


Zi±.\,\ + 1 




Zi 


,2 


Zi 


,2 + 1 







The external face, with h^ ^ 0, is in correspondence of 

open boundaries, so, as we want a single external face, 
if we take periodic boundary conditions in one direction 
(say, along a;), the only possible choice in this framework 
is to take closed and open conditions on the two sides in 
the y direction. 

Our notation is that Zi takes the same value every- 
where: on open boundaries, h~ and 6^ arc determined 
accordingly, while on closed boundaries either we add 
some extra "loop" edges, or we take Zi — t^a > 0. 

The cylindric geometry has all the non-trivial features 
of the original ASM for what concerns group structures, 
polynomial bound on the relaxation time in the group 
action, connection with spanning trees, and so on (even 
the finite-size corrections to the continuum-limit CFT 
are not more severe on a x Ly cylinder than on a 
Lx X Ly open-boundary rectangle) , but, for what concerns 
the determination of Id^ through relaxation of kidf, the 



system behaves as a quasi-unidimensional one, and Mr is 
in general trivially determined. 

Furthermore, in many cases Idr is just the maximally- 
filled configuration Cmax- This fact is easily proven, ei- 
ther by checking that Idr ® Idf = Idr, which is easy in 
the quasi- 1-dimensional formulation (or in other words, 
by exploiting the translation symmetry in one direction) , 
or with a simple burning-test argument, in the cases in 
which on each site there is a single incoming edge from 
sites nearer to the border. Figure 3 shows a few exam- 
ples. 

Such a peculiar property has a desiderable consequence 
on the issue of inversion of a recurrent configuration. In- 
deed, if C and C are in R, wc have that 

(-C) ® C" = Idr ® (-C) © C" = {Idr - C) © C" (10) 

where, if C = {zi}, (— C) = {zi} is the inverse configu- 
ration we seek, but Idr — C is simply the configuration 
with Zi = Zi — Zi — 1, and has all non- negative heights. 

PSEUDO-MANHATTAN LATTICE 

From this section on, we will concentrate on square ge- 
ometries, on the square lattice with edges having a given 
orientation, and all vertices having in- and out-degree 
equal to 2. 

For this reason, we can work with Zi =2, so that S = 
{0,1}^ instead of {0,1,2,3}^ (a kind of simplification, 
as from "CMYK" colour printing to "black and white" ) . 
In this step we lose in general a bit of symmetry: e.g. on 
a square of size 2L we have still the four rotations, but 
we lose the reflections, which are arrow-reversing, while 
on a square of size 2L -|- 1 we lose rotations of an angle 
7r/2, and only have rotation of n. 

Square lattices with oriented edges have already been 
considered in Statistical Mechanics, especially in the 
variant called "Manhattan Lattice" (i.e. with horizontal 
edges oriented east- and west-bound alternately on con- 
secutive rows, and coherently within a row, and similarly 
for vertical edges), cfr. for example [14]. However this 
lattice in two dimensions is not a planar alternating di- 
rected graph, and the results for the related ASM model 
will be discussed only briefly in last section. 

We start by analysing a less common variant, which 
is better behaving for what concerns the ASM model, 
and which we call pseudo-Manhattan lattice (it appears, 
for example, in the totally unrelated paper [15]). In this 
case, the horizontal edges are oriented east- and west- 
bound alternately in both directions (i.e. in a chequer 
design), and similarly for vertical edges. As a result, 
all plaquettes have cyclically oriented edges. A small 
portion of this lattice is shown in Figure 1, where it is 
depicted indeed as the prototype planar alternating di- 
rected graph. Conventionally, in all our examples (unless 
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Figure 3: Reduction to a quasi- 1-dimensional system for the 
ASM on cylindric geometry, on a few examples all having 
Idr coinciding with Cmax- Prom top to bottom: a portion of 
the square lattice, in the two orientations; of the hexagonal 
lattice, in the two orientations, of Manhattan and pseudo- 
Manhattan lattices. Plain edges correspond to Aij = Aji = 
— 1, while a directed edge (from i to j) with k arrows corre- 
spond to Aij = —k. In-(out-)coming bold arrows correspond 
to b~ (6^) equal to 1, while the lozenge-shaped double-arrows 
correspond to 6~ = 1. A loop with k arrows on i corre- 

sponds to Zi — Au = k, otherwise Zi = An. 



otherwise specified) we fix the orientations at the top- 
right corner to be as in the top-right corner of Figure 1. 

Quite recently, in [19] both Manhattan and pseudo- 
Manhattan lattices have been considered as an interest- 
ing variant of the ASM model (the latter under the name 
of F-latMce), with motivations analogous to ours. This 
corroborates our claim that these variants are natural 
simplifications of certain features in the original BTK 
model. 

The main feature of the Creutz identities, on square 
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Figure 4: The top-right corners of the identities Idi'"^ for 
L = 21 and L = 64 (remark: 64 = 3 • 21 -|- 1). The smaller- 
size one is exactly reproduced at the corner of the larger one, 
while the rest of the latter has an evident regular structure, 
according to theorem 1. 



portions of the pseudo-Manhattan lattice with even side, 
is self-similarity for sides best approximating the ratio 
1/3, up to a trivial part, as illustrated in Figure 4. 

The precise statement is in the following Theorem 1. 
Call Id^^^ the set of heights in the Creutz identity for the 
square of side 2L, encoded as a. L x L matrix for one of 
the four portions related by rotation symmetry. Say that 
index (1, 1) is at a corner of the lattice, and index {L, L) 
is in the middle. We have 



Theorem 1 Say L = 3^ -I- s, with s = 1,2,3. Then, 
Id[^^ is determined from Id^^^ and a closed formula, and 
thus, recursively, by a deterministic telescopic procedure 
in at most [logs -^J steps. 

Fors = l or 3, we have {Idi'^^).. = {ld^/^)ij ifij < 
otherwise (id^^^).^ = iff, for i+ j + s even. 



i<£, 

i < 2e, 



\2£- 



j\-i>0; 

<2e- 

j = L- 



(11a) 
(lib) 
(He) 
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and, for i+ j + s odd, 

i>£, 
i>£-l, 



\L-e-i\-j>0- (lid) 
\2£ + s-j\+i<2e; (lie) 



If s = 2 the same holds with i and j transposed in Id[^^ 
and Idf^ (hut not in (11)). Furthermore, Ul = ^^^^^^ . 

The statement of equations (11) is graphically repre- 
sented in Figure 5. 

The understanding of the Creutz identity on square 
portions of the pseudo-Manhattan lattice is completed 
by the following theorem, relating the identity at side 
2L -h 1 to the one at side 2L. We encoded the identity at 
even sides in a LxL matrix , exploiting the rotation 
symmetry, such that the extended 2L x 2L matrix has 
the property 



{Idi%,j = {Id'^'^>)2L+l-J,i. 



(12) 



We can similarly encode the identity at odd sides in a 
structure of almost 1/4 of the volume, namely a {L + 

1) X (L + 1) matrix Id[. , using the fact that 

= {ldf^)2L+2-i,j = {idf^)ifiL+2-j ■ (13) 

Theorem 2 We have {Idf^)i,j = {Id[^\j ifiJ < L, 

{Idl ^)L+i,j =0 if L ~ j is odd and 1 if it is even, and, 

for i < L, {Idj. )i,L+i =0 if L — i is even and 1 if it is 
odd. 



The statement of this theorem is illustrated in Figure 6. 

The proof of Theorem 2, given Theorem 1, is easily 
achieved through arguments completely analogous to the 
ones of [4, sec. 7]. We postpone the (harder) proof of The- 
orem 1 to the discussion of the equivalence with Theorem 
4 below. 

There are some similarities and differences with the 

identities in the BTW model. In that case, besides the 
evident height-2 square in the middle, there arise some 
curvilinear triangolar structures, mostly homogeneous at 
height 3, but some others "texturized" , so that there 
exists a variety of patterns which appear extensively at 
large sizes (in the metaphor of CMYK offset printing, like 
the way in which composite colours are produced!). A 
first attempt of classification of these structures appears 
in ??. In our case, we only have "black and white", but, 
as four zeroes in a square arc a forbidden configuration 
(as well as many others too dense with zeroes), we can 
not have extesive regions of zeroes (black, in our draw- 
ings of fig. 4). A big square in the middle is still there, 
rotated of 45 degrees, while triangoloids are replaced by 
exact "45-45-90" right triangles in a "texturized gray" 
coming from a chequered pattern. Indeed, to our knowl- 
edge, such a regular structure as in theorem 1 was not 



deduciblc a priori (in particular, not before the publica- 
tion of [19]), and our initial motivation was to study the 
emergence of patterns in a 2-colour case. 

The statement of Theorem 1 would suggest to look for 
similar features also in the BTW model. It turns out 
that, while in the directed case the [(L — 1)/3J size is 
fully contained in the corner of the L size, in the BTW 
the \L/2\ size is partially contained in the corner of the 
L size, in an empirical way which strongly fluctuates with 
L, but is in most cases more than 50% (cfr. figure 2 for 
an example). 

The heuristics in the case of the Manhattan lattice, 

under various boundary prescriptions, are somewhat in- 
termediate: the continuum-limit configuration exists and 
coincides with the one for the pseudo-Manhattan, but the 
[(L — 1)/3J corners are only partially reproduced by the 
size-i? identities. 

How could one prove Theorem 1? An elegant algebraic 
property of the identity is that it is the only recurrent 
configuration for which all the "charges" are zero (see 
[4], eqs. (3.3) and (3.4)). More precisely, and in a slightly 
different notation, given a whatever ordering of the sites, 
and any site z, and calling Aij the minor {i, j) of a matrix 
A, we have 



Qj{C) :=^z,(-l)^+^ detAi 



(14) 



and Qj{Idr) = g^det A with qj e Z, with Mr being the 
only stable recurrent configuration with this property. 

A possible proof direction (that we do not follow in 
this paper) could have been as follows. An algebraic re- 
statement of the expression for the charges is achieved in 
Grassmann calculus, through the introduction of a pair 
of anticommuting variables -ipi, ipi per site. Then, by 
Grassmann Gaussian integration, we have that 



(-l)*+^det Aij 



(15) 



where there is a contribution b^ijjkipk in the exponen- 
tial for each site on the border, and a contribution 
{•iph — '4>k)i^k if a particle falls into h after a toppling 
in k, i.e. = — 1 (wc are pedantic on this because the 
asymmetry of A could create confusion on who's who 
with ?/) and tf)). 

So the configurations in the same equivalence class of 
the identity are the only ones such that, for each j, the 
"expectation value" 



(xi^'^O^j) '''' 



JVii;,i;){j:,Zii;^)i;,e^^^ 



JV{ilj,ij)e 



(16) 



is integer-valued. 

Furthermore, expressions as in the right hand side of 

(15) are related, through Kirchhoff theorem, to the com- 
binatorics of a collection of directed spanning trees, all 
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Figure 5: The non-recursive part of the identities Idi^^ (top-right corners) for L = 3£ + s, and s £ {1, 2, 3}, illustrating the set 
described by (11). For s = 2, we show the transposed of /rfr- Black and white stand respectively for Zi = and 1; gray regiorns 
corresponds to chequered parts, starting with white on the cells cutted by 7r/4-inclination lines. 



rooted on the boundary, with the exception of a single 
tree which instead contains both i and j, and the path 
on the tree from i to j is directed consistently. 

So, a possible approach by combinatorial bijections 
could be to prove that, for any j, there is a suitable cor- 
respondence among the forests as above, and a number 
Qj of copies of the original ensemble of rooted forests. 

Such a task is easily performed, even for a generic (ori- 
ented) graph, for what concerns the frame identity Idf, 
for which all charges qj arc 1. Unfortunately, for what 
concerns Idr, and with an eye to the proof performed 
in the following section, it seems difficult to pursue this 
project at least in the case of the pseudo-Manhattan lat- 
tice on a square geometry. Indeed, if the relaxation of 
fci/df to Idr requires tj topplings on site j, it is easy to 
see that qj{Id^) ~ — tj, and from the explicit expres- 
sions for Ul and the values of the topplings (the latter 
are in the following Theorem 4) we see that the values of 
the charges qj are integers of order i^. 



PROOF OF THE THEOREM 

Here we perform the direct proof of Theorem 1. It is 
fully constructive, a bit technical, and maybe not spe- 
cially illuminating for what concerns the algebraic as- 
pects of the problem, but still, it makes the job. 

Now, in order to better exploit the geometry of our 
square lattice, we label a site through a pair ij denot- 
ing its coordinates. We can introduce the matrix t'>^\ 
which tells how many topplings site ij performed in the 
relaxation of k^Idi into Id-^. Exploiting the rotational 
symmetry, we take it simply L x L instead of 2L x 2L, 
with = (1, 1) for the site at the corner, analogously 
to what we have done for {Idi^^)ij. Of course, although 
we call T^^'i and Id[^^ "matrices", they have a single 
site-index, and are indeed vectors, for example, under 



the action of A. 

Clearly, T is a restatement of Idr, as 



kLb± 



^3->^ 3 i' j' 



(17) 



so that Idr is determined from T, but also vice-versa, as 

A is invcrtiblc. Actually, through the locality of A, one 
can avoid matrix inversion if one has some "boundary 
condition" information on T, and the exact expression 
for Idr, e.g. if one knows T on two consecutive rows and 
two consecutive columns (and we have a guess of this 
kind, as we discuss in the following). 

The constraint that {Idr)ij G {0, 1} gives that T 
is locally a parabola with small curvature. More- 
over, in the regions corresponding to homogeneous 
portions of Idr, T must correspond to a discrctized 
parabola through easy formulas. The telescopic nature of 

{Idr > ^'^r ^ , . . .) has its origin in an analogous state- 
ment for (T(^), T(L-^J), . . .), and on the fact that fcz, has 
a simple formula. As we will see, these facts are easier to 
prove. 

We start by defining a variation of T which takes in 
account explicitly both the height- 1 square in the middle 
of Idr, and the spurious effects on the border. Define 
k{L) = [{L — 1){L + 2)/4j (which is approximatively 
khj'i), and introduce 

ff := Tf - - (^-+1) - (^1+1) . (18) 

A first theorem is that 



Theorem 3 



f^^^ = 

^ < 



fori + j>L + (.; 
for all i, j. 



(19a) 
(19b) 
(19c) 
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Figure 6: The recurrent identity on the pseudo-Manhattan 

lattice of side 50 (up) and 51 (down), an example of how the 
identity on side 2L + 1 is trivially deduced from the one on 
side 2L. 



This implicitly restates the claim about the middle 
square of height 1 in Id^, and is in accord with the upper 
bound on the curvature of T given by {Idr)ij < 1- With 



abuse of notations, we denote by t^- 



matrix corresponding to 



also the L x L 

in the £ x i corner with 



i,j < and zero elsewhere. Then the rephrasing of the 
full Theorem 1 is 



Theorem 4 // M^^^ = -(T^ ^_ . - T^^^, ^_ .) for s = 
1,3 and the transpose of the latter for s = 2, defin- 
ing 0{n) = 1 for n > and otherwise, and the 
"quadratic + parity" function q{n) on positive integers 



g(n) = ein) (n^ + 2n+ L_LDl^ 



n < 0: 



(20) 



(n + 1)^/4 n positive odd; 
n{n + 2)/4: n positive even; 



dratic + parity" on a finite number of triagular patches 
Mlf^ =e{i-L + l)^- (i) - max(0,i -l-l) 
+ q{j -i + L-£)-q{j + i-L-e-4) 

+ e{j-L + e)[-{i)-e{L-e-i + i) (21) 

+ qii-j + L-e-l) + qii+j-L + e-2) 
+ q{i+j-L-£-3) 

Clearly Theorem 3 is contained in Theorem 4, just by 
direct inspection of the summands in (21). The equiva- 
lence among Theorems 1 and 4 is achieved still by direct 
inspection of (21), with the help of some simple lem- 
mas. First remark that the use of — T instead of T 
makes us work "in false colours", i.e. effectively inter- 
changes z into 1 — in Mr- Then, defining Vxf{i,j) ■= 
/ (« + 1 , J ) + / - 1 > j ) - 2/ (« , J ) > and analogously with 
±1 on j, we have that 



V',^a(±*±J-a) = 



1 ±i ± j — a > and even; 
otherwise; 



(22) 



then M(^) is a deterministic function, 



piecewise qua- 



and that V^max(0,j — a) = 1 at j = a only, while 
Wy max(0, j — a) = always (this reproduces the £ extra 
zeroes out of the triangles depicted in figure 5). 

So, the explicit checks above can lead to the conclusion 
that equation (17) holds at every L for the expressions 
for T and Idr given, in a special form, i.e. subtracting 
the equation for size L with the one for size £, but on the 
Lx L system. The latter is then easily related to the one 
on the £x £ system, at the light of equation (19a), which 
allows to state that there are no different contributions 
from having splitted the four quadrants of T^^\ Thus, 
equation (19a), guessing the exact expression for Tij on 
the edges of each quadrant by mean of a simple formula, 
is crucial to the possibility of having the telescopic recon- 
struction procedure. 

In other words, even without doing the tedious checks, 
we have seen how, by a series of manipulations corre- 
sponding to the subtraction of £x £ corners to the Lx L 
matrices, the proof is restricted to the analysis of the 
"deterministic" part, with i > £ or j > £. As, in this 
case, all the involved functions depend on L through its 
congruence modulo 2, or 3, or 4, it suffices to check nu- 
merically the theorems for 12 consecutive sizes in order 
to have that the theorem must hold for all sizes. We did 
the explicit check for sizes up to i = 64. 

The only thing that we need in order to conclude that 
the conjectured expression Id^ corresponds to the iden- 
tity is to prove that it is indeed a recurrent configura- 
tion. Again, we do that in two steps, in order to divide 
the behaviour on the self-similar £ x £ corner and on the 
deterministic part. 

First, remark that in the internal border of the 

quadrant (the sites ij with i ov j = L) is burnt in the 
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Burning Test even without exploiting the other mirror 
images. This is true because, at every L and on both 
borders, we have a substructure of the form 



which is burnt through the foUowing cascade (we put the 
burning times, and denote with arrows from a; to y a 
toppUng occurring in x which triggles the toppUng in y) 



2 ;.3 



5 -: 4 5 



6 >7 < 6 >7 



9 < 8 >9 -: 8 >9 12 H3 



10>11-:10>11.:10>11 14:.15 18:49 22:--- 



13<12:.13<12:.13<12il3 16:47 20:-21 



We choose to show just an example, instead of giving the 
explicit formula, but it should be clear from the pictures 
above that a gcneral-L regular procedure exists. 

So, Mr satisfies the burning test if and only if the de- 
terministic part of Idr (the L x L square minus the £x £ 
corner) satisfies the burning test with the border of the 
corner having a bij = 1 every two sites. This check must 
be performed only on the deterministic part, and again is 
done in a straightforward way, or, with conceptual econ- 
omy, implied by the explicit numerical check on 6 con- 
secutive sizes. This completes the proof of all the three 
theorems. □ 



MANHATTAN LATTICE 

For the Manhattan lattice, we performed an investiga- 
tion similar to the one above for the pseudo-Manhattan 
one, although, as we already said, the motivations arc 
less strong, as this lattice is not planar alternating. A 
small-size example of Creutz identity, compared to the 
pseudo-Manhattan one, is shown in Figure 7. 

The numerics gave positive and negative results. The 
positive results concern the filling number k^, that, ac- 
cording to extensive tests (up to i ~ 100) we conjecture 
to be 



kL = 



_] \L{L + 2) 
■{L + lf 



L even 
L odd 



(23) 



Furthermore, the whole quadrant except for the Ix (€— 1) 
corner seems to be deterministically described by a set 



■11 
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Figure 7: The recurrent identity on a portion of side 50 of the 
square lattice, with pseudo-Manhattan (up) and Manhattan 
(down) orientation. 



of rules analogous to equations (11) (again, £ = \_{L + 
1)/3J), and, as in (11), depending from the congruence 
of L modulo 3, and a transposition involved if s = 2. 

The negative result is that the telescopic exact self- 
similarity between side L and side £ seems to be lost in 
this case. As we said, this property relies crucially on the 
simplicity of the toppling matrix on the boundary of the 
quadrants, which seems to be an accidental fact of the 
pseudo-Manhattan lattice, and has few chances to show 
any universality. For this reason we did not attempt to 
state and prove any theorem in the fashion of Theorem 1 
in this case. 



CONCLUSIONS 

We have studied, numerically and analytically, the 

shape of the Creutz identity sandpile configurations, for 
variants of the ASM, with directed edges on a square 
lattice (pseudo-Manhattan and Manhattan), and square 
geometry. An original motivation for this study was the 
fact that heights are valued in {0, 1} (while the original 



BTW sandpilc has heights in {0, 1, 2, 3}), and wc conjec- 
tured that this could led to simplifications. The results 
have been even simpler than expected, and qualitatively 
different from the ones in the BTW model. 

In the BTW model, the exact configuration seems to be 
unpredictable: although some general "coarse-grained" 
triangoloid shapes seem to have a definite large-volume 
limit, similar in the square geometry and in the one ro- 
tated by 7r/4, here and there perturbations arise in the 
configuration, along lines and of a width of order 1 in 
lattice spacing. We discuss the role of these structures in 
various aspects of the ASM in a forthcoming paper [13]. 

The triangoloids have precise shapes depending from 
their position in the geometry, and are smaller and 
smaller towards the corners. Understanding analyti- 
cally at least the limit shape (i.e. neglecting all the sub- 
extensive perturbations of the regular-pattern regions) is 
a task, at our knowledge, still not completed, although 
some first important results have been obtained in [18], 
and further achievements in this direction have come with 
the work of Levine and Peres [16, 17], both in the similar 
context of understanding the relaxation of a large pile in 
a single site (see in particular the image at page 10 of 
Levine thesis, and the one at 
http : //math . berkeley . edu/ 

-levine/ gallery/ invertedscindpilelmlOx . png, 
and, for the directed model, the one at page 20 of 
newt on . kias . re . kr/ 

-nspcsOS/Presentation/Dhar .pdf ). 

In our Manhattan-like lattices on square geometry, 
however, we show how the sitiiation is much simpler, and 
drastically different. Triangoloids are replaced by exact 
triangles, all of the same shape (namely, shaped as half- 
squares), and with straight sides. All the sides of the 
triangles are a fraction ^3"*^ of the side of the lattice (in 
the limit), where the integer fc is a "generation" index de- 
pending on how near to a corner we are, and indeed each 
quadrant of the configuration is self-similar under scal- 
ing of a factor 1/3. The corresponding "infinite- volume" 
limit configuration is depicted in Figure 8. A restatement 
of the self-similarity structure, in a language resembling 
the z —>■ conformal transformation in Ostojic [18] 

and Levine and Peres [16], is the fact that, under the 
map 2 — > Ina;, a quadrant of the identity (centered at the 
corner) is mapped in a quasi -doubly-periodic structure. 

Also the filling numbers (i.e. the minimal number of 
frame identities relaxing to the recurrent identity) have 
simple parabolic formulas, while in the original BTW 
model a parabola is not exact, but only a good fitting 
formula. 

These features reach the extreme conseqiicnccs in the 
pseudo-Manhattan lattice, where the exact configura- 
tion at some size is deterministically obtained, through 
a ratio- 1/3 telescopic formula. These facts are not only 
shown numerically, but also proven directly in a combi- 
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Figure 8: Up: the limit Creutz Identity configuration on our 
Manhattan and pseudo-Manhattan lattices. White regions 

correspond to have height 1 almost everywhere. Gray regions 
correspond to have height and 1 in a chequered pattern 
almost everywhere. Down: the image of a quadrant under 
the map a — » In 0. 



natorial way. 
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